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1.  Introduction. 

In  a  number  of  applications  It  la  necesaary  to  knov  the  propertlea  of 
the  product  of  random  variables;  this  occurs  In  particular  when  the  random 
variables  Involved  have  dimensions  of  a  ratio  like  tolerance  expressed  In 
percentages  of  desired  values,  fuel  consumption  per  mile,  amplification 
ratios,  etc.  The  special  situation  discussed  in  this  paper  of  the  product 
of  a  number  of  Independent  Identically  distributed  random  variables  arises 
for  Instance  In  the  case  vhen  some  devices  designed  to  amplify  a  magnitude 
and  having  identical  characteristics  are  connected  In  series.  If  art>  Is  the 
random  variable  describing  the  amplification  by  the  i  th  device  the  total 
amplification  Is  also  a  random  variable  and  It  la  Important 

to  knov  the  distribution  of  this  product;.  Examples  of  engineering  applications 
involving  products  and  quotients  of  random  variables  can  be  found, for  Instance, 
In  Donahue  [4], 
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Zt  vu  shown  by  Springer  and  Thompson  {9]  bow  to  obtain  tha  probability 
danslty  function  (p.d.f. )  of  product!  of  iv  independent,  Identically 
distributed  randan  variables  by  the  application  of  tha  Mellln  transform; 
they  have  obtained  attng  others  formulae  (in  the  form  of  rapidly  convergent 
Infinite  series)  for  the  p.d.f*  of  the  product  of  k  Independent  normal 
variates  and  y\,  Cauchy  variates  and  treated  aleo  a  special  caae  of  Beta 
variates  [see  below  formula  (7)].  Their  method  Is  a  generalisation  to 
factors  of  a  method  presented  by  Epstein  (5)  for  .  Following  Epstein 
many  authors  applied  the  Mellln  transform  to  the  study  of  distribution  of 
products  and  quotients  of  randan  variables;  a  detailed  bibliography  can  be 
found  In  Springer  and  Thompson  [9]  and  [8];  cf.  also  Kotlarskl  [6]  and 
Zolotarev  [10]. 

The  Mellln  transform  of  a  function  where  x.70  Is  defined  as 

m 

(1)  -  Jx'-'jivix  . 

0 

Under  certain  regularity  conditions  [c.f.  Cour  ant -Hilbert  [3]  pp.  10  3 -10  5] 
this  transform  considered  as  a  function  of  complex  variable  S  admits  an 
Inversion  Integral: 


P«*e  3 


c 

(2)  /w  .  lir / 

c-i*0 

where  the  path  of  Integration  le  a  line  parallel  to  the  Imaginary  axle 
and  to  the  right  of  the  origin. 

An  immediate  consequence  of  the  definition  of  the  Mellln  transform  is 
the  following  multiplicative  property:  If  x  and  y  are  Independent  random 
variables  with  the  p.d.f. 's  fix)  and  $(>j  respectively  and  if 
is  the  p.d.f.  of  2*  x.y  then 

(3)  ]  =  tl Lfi*)h]  . 

Thus  to  find  the  p.d.f.  of  y*  xi  x\  where  xt-  are  Independent  Identically 
distributed  random  variables  with  p.d.f.  fix)  it  is  sufficient  to  find  the 
Mellin  transform  of  fl *)  ,  to  take  its  >tth  power  and  to  find  the  Inverse 
with  the  aid  of  formula  (2). 

The  following  property  of  the  Mellin  transform  will  be  needed  later 
in  the  discussion  of  the  probability  distribution  function  of  the  product: 

If  fix)  is  the  p.d.f.  of  the  random  variable  x  which  has  finite 

second  moment  and  if  — 

F(*J  -  f  f(4r)dit  ^Lx)  »  / f'.tjett  «  i~F(x) 

0  '  A 


J*!  WftWp  £*  «*#  *■»  * 
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then  the  Nellln  transform  of  T(x)  la  equal  to 

w  M[fi~>b]  *  /V; r"i*  -  s-x'ttofirfx'/M*  '  t"n[fivl‘«], 

the  vanishing  of  the  ten  S  xr  f(xj  when  <-t"  following  from  the 
existence  of  the  second  ament. 

It  has  been  shown  by  the  author  of  this  report  [7]  that  the  results 
obtained  by  Springer  and  Thoopson  In  [9]  for  normal  random  variables  can 
be  presented  In  a  somewhat  simpler  fon  and  that  a  general  formula  for 
any  number  of  factors  can  be  given  (still  rather  unwieldy  In  the  case  of 
large  values  of  tv).  It  has  been  also  shown  that  by  a  direct  application 
of  the  Mellln  transform  similar  infinite-series  expansions  can  be  derived 
for  the  corresponding  probsblllty  distribution  functions;  this  is  useful 
since  It  Is  not  always  easy  to  obtain  them  by  the  straight-forward  Integration 
of  the  relevant  infinite  series  representing  probability  densities.  Finally* 
attention  has  been  drawn  to  the  fact  Implicit  in  the  Springer -Thompson 
treatment  that  it  Is  sufficient  to  evaluate  the  formulae  for  the  p.d.f. 
of  the  product  of  Independent  exponentially  distributed  random  variables 
and  from  such  tables  the  p.d.l.  for  products  of  gamma,  normal  and  Velbull 
random  variables  can  be  quickly  evaluated  by  simple  transformations. 
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The  practical  uaafulnaaa  of  the  rcaulta  described  abort  la  Halted  by 
tha  fact  that  all  the  corresponding  distributions  hare  Infinite  rancea 
while  the  phyaloal  devices  to  whloh  our  aathaaatlcal  node  la  are  to  be 
applied  often  hare  finite  charaeterletlca.  Thus,  the  problaa  arise  a  to 
what  extent  the  aethoda  applied  In  the  abort -quoted  papera  can  be  Modified 
to  obtain  alallar  re  suit  a  for  randoa  rarlablea  haring  finite  ranges  and 
the  first  family  of  randoa  variables  which  suggests  itself  and  setaa  to  be 
tractable  In  this  way  Is  the  family  of  Beta  distributions  with  parameters  f>, 
l.e.  the  randoa  variables  haring  p.d.f. : 

(5)  P  C*J  M)  s  Si’*'*)  P*1*  1*/  j 

where  la  the  Euler  integral  of  the  first  kind: 

(6)  -  {x'-'C'-xj1"**  -  ro:n$j/r(f*u. 

notice  that  In  (5)  >>oth  />  and  <y  are  assucaed  to  be  not  less  than  unity; 
for  fi  or  ^  leas  than  1  the  p.d.f.  would  exhibit  an  Infinite  Jump  at 
jr=  o  or  and  the  Inversion  formula  (2)  would  not  be  applicable. 

This  Is  not  a  serious  limitation  since  those  "U”  and  "J"  -shaped  distributions 
are  of  negligible  practical  Importance. 
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Springer  and  Thompson  [9]  gave  an  explicit  formula  for  the  product 
of  tv  Independent  Beta  randon  variables  of  a  special  class  l.e.  of  Beta 
variables  with  parameters  />  ■at/j  ■  1  defined  by  the  p.d.f. 

(7)  <Wji j  **■(*+•)  Oixi/Jot>o, 

vhlch  they  called  "monomial"  distributions;  the  rectangular  distribution 
is  a  member  of  this  class  for  &  ■  0. 

In  this  particular  case  the  application  of  the  Mellln  transform  gives 
an  Immediate  answer.  Thus 

,a\  M  '  UK*4'*  ,)li]  ~  z 

\o)  L 

The  Mellln  transform  of  the  p.d.f.  of  the  product  of  *1  Independent  random 
variables  of  that  kind  is  equal  to 

(9)  H  [  ’  C 

Put  the  inverse  Mellln  transform  of  (  4  *  *  )  K  can  be  found  In  tables 
[see  e.g.  Bateman  Manuscript  [2],  formula  7.1  (16) ]  so  that 

,  L  -  log  xj*  '  L  ^(KjJ  If  '<x£tJ 

(10)  (JH 

[  0  it 


1  . 
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There  ia  a  wlde-apread  belief  among  working  statisticians  that  the 
problem  of  the  distribution  of  products  of  Independent  random  variables 
Is  not  of  a  very  great  interest  since  "one  can  always  find  the  distribution 
of  y  ■  log  x ,  apply  the  theorems  on  the  addition  of  Independent  random 
variables  yt-  ■  log  and  then  revert  to  the  product  of  the  x's. "  It 
happens  so  that  In  the  above  particular  case  this  procedure  works  and  the 
application  of  the  Mellin  transform  (although  very  attractive)  Is  not 
essential.  Indeed,  If  we  put  >■  log x  then  it  follows  from  (7)  that 
the  distribution  of  y  given  by  the  p.d.f.  fiyj  is  the  negative  exponential 
distribution  with  parameter  A 

—  yCK+iJ 

j-iyj  -  £ 

and  It  is  well  known  that  the  sum  of  independent  random  variables  of 
this  kind  is  a  gamma  distribution  vlth  the  shape  parameter  K  and  the 
scale  parameter  °C+/  so  that 

/K  „j = (**v  ,  •</*-. 

Putting  x.  ■  exp  (-/)  we  obtain 

k  t / ./  -  <* ♦  o' ' x< i • - ■ l<* *  j"  'lr‘‘ KJJ "j 

which  agrees  with  (10).  However,  this  seems  to  be  an  exceptionally  simple 
situation  and  in  more  general  cases  the  passage  through  the  distribution  of 
the  sum  of  logarithms  does  not  make  the  problem  any  easier. 
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2.  Probability  density  function  of  the  product  of  Beta  variables. 


2.1  Simplification  of  the  problem 

The  Me  11  in  transform  of  (5)  la  equal  to 


(11) 


a  r(w-v  n.b+f)/Ll'iwtt~tJ- 


Denoting  by  the  p.d.f.  of  the  product  of  n.  such  Independent 

randan  variables  ve  have,  in  view  of  (3): 


(12) 


■]  -  /•  W'W ,  lr> rW 


and  (2)  yield.:  J—  Ut  __  . 

r'<n>  *" 

On  substitution  t  *  i+»-/  this  becomes: 

r\nt4)  *H'‘j /  t'>0^ 

(13)  **•  C -im 


the  path  of  Integration  being  still  a  line  parallel  to  the  imaginary  axis 
and  to  the  right  of  it  since  /■>  >  1.  But  the  above  formula  can  be  rewritten 

<  *  '  C 


as 


r _  x/*-/  r2ii*±jj  ,  00 

r^r^nj  2jri  c-/«* 


which,  in  view  of  (13)  with  />  ■  1  gives 

(1U)  )  -  ///v'Y*">, 


f 
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This  shovs  that  it  Is  sufficient  to  evaluate  the  p.d.f.  of  the  product  of 
x.  randan  variables  of  a  family  of  one -parameter  Beta  variables  having 
/»  -  1  and  arbitrary  f,  1  j  for  p  >  1  the  relevant  p.d.f,  are  obtained 
by  the  above  simple  transformation.  In  the  subsequent  tvo  subsections  we 
shall  discuss  the  derivation  of  the  p.d.f,  / $+.1 separately  for  the 
case  when  *y  Is  an  Integer  and  the  case  when  it  Is  not  an  Integer. 


2.2  The  case  of  an  Integer 
From  (12)  ve  have 


(15)  " LI 


»J/ 1  j  =  r'ilJ  /r’'nfSJ  ■ 


But  with  an  Integer  <j  , 
rL<j+tj  *  + 


so  that 


1f  i  ix  -i  >j  i  *  Ls*k’'/ 

\  >  J  j  4  J  k*o  J 


and  the  inverse  formula  yields 


■  .  rr 


* tlLS4icJ 


L  > 0 


i 
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The  Integrand  In  (17)  haa  <y  poles  of  the  Hth  order  for  S  -  -j 
(j  *  Ojtj  . .  j  f-/  j  and  the  Integral  can  be  evaluated  by  contour  integration 


yielding 


(18)  jKl*'- 


r*ti*o 


for  0  ^  X  / 


for  X}/  I  j 


where 


(19)  %( 


(p-0!  (  f  7  /  <'7 


is  the  residue  of  the  Integrand  of  (17)  at  the  Hth -order  pole  at  f  -  -  ’  . 
[For  the  details  of  contour  integration  see  Appx/para  (a),  (b),  (c)J. 

Let  us  write  r, i  for  the  product  of  (w-*/)  factors  -77  *i(Kj 
k**  A'* 

(i.e  the  product  in  which  the  y  th  factor  has  been  omitted)  and  similarly 
V(y/ 

w  .  for  the  sum  in  which  the  /  th  term  has  been  omitted.  Denoting  by 

£»•  ' 

GjiiJthc  function  in  the  brackets  in  (19)  we  have 

T-\  . 


V  i  S* 
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Before  attempting  to  give  a  formula  for  arbitrary  integer  7  let  ua 
discuss  two  simple  cases  of  f  ■  1  and  j  ■  2.  For  7  ■  1  ve  have  to  consider 
only  one  pole  and  Its  residue.  Here  (20)  reduces  to  * ,  and  (19) 


gives 


%lXj  *(%•//.'  *  V“l0«xJ  -  <**/,/  /-iogjcj  / 


so  that 


(21)  /K  ix*  1‘  LJ  " 


'  106  X.J 

'K-ir  v.  ' 


for  0<*Ld-t 


for  l  . 


From  (14) 


(22)  ;UL*‘r>'1)  • 


( r lor  ■LJ 


for  ? 


for  A**-, 


which  agrees  with  (10)  /hen  ft 

For  ^  ■  2  we  have  to  consider  two  residues  for  y  »  0  and  j  ■  1  at  the 
poles  J  >  0  and  4&-1.  Formula  (20)  yields 

,  £,c/' 

while  formula  ( 19)  gives 

A..*/  0)  =  tr7y'  3^-  / <‘V* "-■*//  J  . 

-i.  . 

A.  (x;  n,  IJ  m  ‘  /«■-'. 

Applying  the  Leibniz  formula  for  the  (**>-•  )th  derivative  of  the  product  we  find 

r— 7  *-/  ,  A  ~lK+k+/ 

/  ^  _lag  *  V//'7 
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and 


Hence 


A(k «,  H  *  t*-v'  Z.c**>"r^  i{~0-,*-*** 

V'  ,-,  -Urtx/t-/*' 

'  i  /  I 


l 


(23) 


41*/'  *< 


c*0 


for  x^/  . 


Thus,  In  the  two  cases  ^  ■  1  and  y  ■  2,  an  explicit  formula  for 
i  vi4;  4  yj  has  been  given.  For  higher  values  of  ^  the  application  of  the 
Leibniz  rule  would  lead  to  formulae  which  would  become  more  and  more 
complicated  and  It  is  suggested  to  take  advantage  of  the  fact  that  the 
logarithmic  derivatives  of  (20)  can  be  easily  obtained;  once  they  are  found 
the  problem  is  solved  since  it  is  known  how  to  obtain  the  >vth  derivative 
of  a  function  from  its  logarithmic  derivatives.  Indeed,  if 


(24)  -  -  [log  ^ 


and 


*r>  JL  >/; 

i *  dsr  *l  J 


then 


(25) 


i  ’"t-v 


't‘j .  ZK  [ti'i.  .  ■  .,&*>], 


where 


j  r’j 


(26) 


\ V X<-*  J  ‘  ■ 


I 
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here  the  itn  li  extended  to  all  the  partitions  of  number  n.  euch  that 

JP)  >.  jtr) 

n.^*»tWvl  denotes  Mi)  and  /f  the  ^derivative 
of  /fy'Jvith  respect  to  i  . 


Clearly  , 

(26a) 

Z,  =  4- . 

(26b) 

z  = 

^*2-  ; 

(26c) 

-  /f  V  -f 

4 

and  further  formulae  can  be  derived  recursively  by  applying 

(„)  ZK+,  = 

(Cf.  for  example  [7]). 

In  our  case  is  given  by  (20)  and 


(28) 


'or,  f,  (*)  '  -  S  '■  3  *  "  >v  21 
7  A*o 


lot'  ( 
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The  logarithmic  derivative  of  (20)  is 

Ajii)  *  -  log  r  - 
*  *  •* 


and 


-(rH) 


/V  J 


Since,  according  to  (19)  theae  derivatives  are  to  be  taken  for  X  -  — 

#-/ 

ve  have  V-i,/  ,  _ *, 

^  -  *'  £  ,  4_, 


(29) 


Ajrij  "  ~  log  •*■  '  *  TZJ(-k'jJ  - 


and  the  final  formula  for  the  p.d.f.  of  the  product  1b  ‘ 


(30)^ 


-  Fa'  '  y^Uy/'Z-,  [ \  4'...,  4**], 


where  the  arguments  of  the  Z  -function  are  given  by  (29). 


Example.  In  the  case  of  a  rectangular  distribution  ^  ■  1, 
los  .  -j  Iok  *  ^  «  -1°6  Jt, 

all  the  derivatives  of  4t  l  V  with  respect  to  ■>  vani  shing.  Hence  4-, '  A 

and  ^  ^stMA'(-l°gJC),  which  agrees  with  (2l). 

Example.  In  the  case  />  •  1,  ^  ■  2  it  is  easy  to  verify  that 
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4(<l-  -io,x  •iW',  4'l(0’*r'M‘-On' 

S, (t)  »  x‘i'K  ,  A,/iJ  ■  -  log*  .ns"  ,  A*(l)rn 

Substituting  £ ■  0  lo  the  first  line,  £  ■  -1  in  the  second  we  have 

-  i  ,  A.l°)  - 10«  Jt  -  x  ,  -  *  rV'T" 

C,H  «  xHJK,  A,M9  ‘  lo<?*  *  4?^*. )  *  nr/ 

Pomula  (90)  and  (26a)  gives  for  »v«  2: 

(31)  /,y  =  **7-  >**  -  V  **f***  +  ij} , 

and  fonsula  (90)  and  (26b)  gives  for  x  ■  3: 

(32)  ft j(JLj  *.*-J  =  -^/(‘logjc  +  if’-hi]j  ; 

these  results  can  be  verified  by  putting  “H»  ”,  X  *  ,j  etc.  in  (2J). 

2.3  The  case  when  g  is  not  an  Integer. 

Fran  (12) 

(33)  i‘Pli] 

end 

(3U) 

Here  the  integrand  has  an  infinity  of  poles  of  the  \th-order  at  S  ■-/ 

(  *  *  0,1,...)  and  the  Integral  in  (3^)  can  be  evaluated  by  contour 
Integration  yielding 


r\i+o 


i-.  f 

>.7u  J 


t  *>\j» 


r(inj 


K 

< 


Ui 


C "?  0 


for 
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(35) 


"iro  jL  KL 


*jjj 


for  X  >/  L  t 


where 


,  un"  I  ,-j  CliiL  u*j\Kl 
=  (py  Zr-  l  r'h*‘]  J  J  ‘‘-j 


la  the  residue  of  the  Integrand  of  (34)  at  the  nth-order  pole  at  /  •  -j  . 
(^or  the  details  of  contour  integration  see  Appx/ paras  (d)  -  (j) ] .  Clearly 

n*n**h~  1  ri.*w£tw 

ami  for  )  ■  0  the  above  Identity  is  still  valid  if  we  agree  to  interpret 
the  "empty"  product  as  unity.  Hence 


(37)  Rl*‘ 


:i  *£l'lx 'srnO^*tjr  (stkJ  j  / 

*'KJJ  *  IP  C  i%°  J /** 


V  . 


Denoting  by  £jts)  the  function  in  the  brackets  in  (37) 
we  have  .  ilL  .  _  k 

Cj[s)  *  x‘r"(‘*j*'J  r  , 

log  ■  -  1  log jc  -#■  n  log -  Klogf(VfJ*nz^” log  [s+k  J 


— W 

(interpreting  again  the  "empty"  sum  2^  as  zero  to  retain  the  validity 
of  the  above  for  J  ■  0. ) 


The  successive  logarithmic  derivatives  of  are 


(38) 
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Aj(o «  -  log  x  *  k ~ *■  td - K± ("W1 
J  «  %?>(>*)*>) +t-d”r.' i- 

where  yi)  denotes  the  logarithmic  derivative  of  ?*(•)  (the  Euler  "P8i,,-funotlonJ 
lr\. 

and  Y'  c  J  its  eueceeslve  derivatives. 

Consequently  J-  /  XJ(r,V 

GjH)  c  rlliJ  ~ (j!j'rKhy)  J 

09)  Ajl'il  =  -  108  *  *  *)**>  "•  rir/J  *  n  h:^" J 

and  by  the  argument  which  has  been  applied  in  derivation  of  fonnula  (30) 


we  have 


(1*0)  PkLxS  * 


where  the  values  to  be  put  into  ^./  -function  are  given  by  (39)* 


The  infinite  series  (40)  is  absolutely  convergent  for  0  <*<  1.  Indeed, 

the  coefficient  ni«)0V  rt  1-j)  is  equal  to  the  K  th  power  of 
1C1-V  ■  to-j) 


ft/  m  Y\  . .  rj*  (^“0,1 ,...).  For  y  -  0,1,  there  are  [}]  «■  1 

* 

positive  values  of  ;  let  ^  be  their  maximum.  For  J  >  /f  7+  j  ^  j~ 

f  <1  80  that  /  <  <?  *  Further,  the  values  to  be  put  into  the 
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„  ,  -function  are  given  by  (39)  and  since  ,  i«  not  an  integer  the  aeries 

expansions  of  *p°Cy)  valid  for  yi°, .  can  be  applied  to  evaluate 

[see  e.g.  [l]  formula  (6.4.10)).  Thus  /  VOf-Jj/  -  r  *  r.'fo-tilj  J. 

This  is  bounded  from  above  by  a  constant  independent  of  j  and  the  sane 

applies  to  the  remaining  terms  of  (see  e.g.  [l]  formula  (6*4.2)]. 

The  terms  of  Aj  are  not  bounded,  but,  by  a  similar  argument  it  can  be  shovn 
that  they  are  of  the  order  of  0(logy  )  for  any  fixed  K  >  0  [aee  e.g.  [l] 
formula  (6.3.2)  and  the  recurrence  formula  (6.3.6)].  Since  the  values 

(39)  enter  into  the  2h.,  -function  as  the  nth  powers  at  the  most  the  series 

(40)  is  dominated  by  J  multiplied  by  a  constant  and  thl*  still 

shows  that  it  is  convergent  for  positive  X  smaller  than  1. 


3.  Probability  distribution  function  of  the  product  of  Beta  variables. 
General  Remarks. 

3.1  Let  the  probability  distribution  function  corresponding  to  the  p.d.f. 


(5)  be 

(41) 


x 


I  3  (4j  />,f) 


In  view  of  (4)  and  (ll)  we  have 

<«)  V  l1  -  Snt*-'  fi-t'l  ‘J  *  S~'r'(p+!J  , 
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and  the  inversion  formula  (2)  yields 


c#/*» 

(1*3)  .  c,0. 


It  is  not  possible  to  simplify  the  problem  by  reducing  the  discussion 
to  the  one-parameter  family  of  probability  distribution  functions  in  the 
way  similar  to  that  carried  out  in  section  2.1  for  the  p.d.f.  's.  We  have 
to  investigate  the  tvo-parameter  family  of  probability  distribution  functions 
(U3)  and  ve  oh <411  again  discuss  separately  the  cases  of  Integer  and  non-integer 


3.2  The  case  of  Integer  . 

Since,  for  Integer  $ 

FCfi+l  +  Sj  *  ’  •  (S'*1-*  F [p+Uj 

(43)  becomes 


'^u.j \ ,  '->0 

Lb  J  *.»• 


The  Integrand  has  one  single  pole  at  J  ■  0  and  ^  poles  of  the  *1  th  order 
at  f  •-j-h  (/■  0,1,  ).  The  integral  can  be  evaluated  by  contour 

integration  by  the  method  similar  to  that  applied  in  the  evaluation  of 
integral  (17)  [cf.  Appx  para  k].  Since  the  residue  at  J»  0  multiplied  by 
is  clearly  equal  to  1  ve  have 


i 


I 


(1»5)  &K 


-  i\ ■■  v/ 
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for  0<x<  I  ,, 


for  4T>  /  t 


q*  , 

where  a  6<>' V y  It  the  residue  of  the  integral  appearing  In  (44)  at  the  >i 

order  pole  at  f  ■  -/-/•  ( j  ■  0,1,...,^*/): 

<<*)  k*L*.'Kj)  '  cf7T/ "OVrfiJ  /  . 


Retaining  the  notation  Introduced  in  formula  (20)  we  can  write  the 
function  Gj  ^An  the  brackets  of  (46)  as 

f-‘  . 

(47)  0^(1)*  X  S  //  Lf+l+k) 

J  k*o 


Starting  again  with  the  discussion  of  two  simple  cases  we  have  for  f,  ■  1 
only  one  residue  at  J  ■  0.  From  (47)  ■  *  fi_/  and 

)  9  ZiCrJ  *  (-'J  in-r-/j/s~*  J/ 

y  S*-f» 


SO  that  -  .  1.  IT  -nrr 

°cx)  P  ; 

M~'  br  r 

(48)  K'-r  <fr,  (  -logjt.y  . 


r  ~n.tr 


The  derivative  of  the  above  Is  equal  to 

Ik  (*;/>,  iy/  -*  ^  ’og  -c/  /£*  //  y 

which  agrees  with  (22), 


I 
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For  p  ■  2  we  have  to  consider  two  residues  for  /  ■  0  and  j  •  1; 

fC(x;%t0)  -  {xT's  'ifi+iHj'*}  /  [%  m/§  ^ 

0 * (jhj!  itF-  {x~Sr‘ /t*-*-,  . 

Applying  the  Leibniz  Formula  for  the  (  *i~/  )th  derivative  vo  obtain 

--  itry  2LCr)£.ll)  4  -l0«  -  (H  A 


Evaluating  the  above  at  1  ■-/”  and  i  ■  -/>-/  respectively 

a*i<, «,  n  -  -  £<"4  */ 


so  that 


t>  (”lo^v  [x?P  [r*»j  1 


The  derivative  of  (49)  is  equal  to 

«»  / 

/  k.  v.  h  T  1  rT  p  /»-/  r+'  i 

c*o 

and  Is  equal  to  /i t Jj 
in  view  of  (23)  and  (14). 

Again,  ae  In  section  2,2,  for  higher  values  of  the  application  of 
the  Leibniz  rule  would  lead  to  complicated  expressions  and  we  shall  again 
evaluate  the  logarithmic  derivative  of  (47).  Here 


X 
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r*  /  , 

log  ry  (*J  -  _  SlOgX  -  logs  -  *1  log  + 


and  the  logarithmic  derivative  of  (47)  is 


(50)  •  -  l0*x 


while  the  higher  derivatives  are 


/*>  -  r',w/V'"<V  ■ 

(51)  4  <■*>  -  L  J  k-° 


Since  according  to  (46)  these  derivatives  are  to  be  taken  for  J'=7'/,> 


we  have 


C]<ri-N-  £ "W;: 

/.  ...........  a  . 


(52)  Aj  r/-'V  -  --O'*  ->■  jvp  *  J  J  / 

*(+>  .  /  /  ’  •(***)  _u  ml*'  ’  ^  >  ‘  -  l*”*1) 

Aj  L’/~/>)  •-  r-  U*tu  +  £r; 

and  the  final  formula  for  the  probability  distribution  function  of  the 


product  is  given  by 


(53) 


(*;/•/ r; * L*-*/  j  ^*w J/  j  J  ^ ' 


where  the  arguments  of  -  xv  *  function  are  given  by  (52). 
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Example.  If  ^  ■  1,  according  to  (53)# 


Th's  formula  la,  for  high  values  of  k  ,  more  complicated  than  (48)  but 
It  can  be  verified  that  for  n,m  2,3  etc..  It  le  equivalent  to  (1*6). 

3. 3  The  caae  when  la  not  an  Integer. 

Since  ^  la  not  an  Integer  the  Integral  (1*3)  cannot  be  written  in 
a  simpler  form  aa  in  the  caae  of  (1*4).  The  integrand  of  (43)  has  a  single 
pole  at  S»  0,  and  an  infinity  of  poles  of  the  -nth  order  at  1* 

(  j  ■  0,1,...).  The  Integral  in  (43)  can  be  evaluated  by  contour  Integration 
by  a  method  similar  to  that  applied  in  the  evaluation  of  the  integral  of  (34). 
[cf.  Appx,  para  (k)).  Since  the  residue  at  j  -  0  multiplied  by  ^ 


equals  to  1  we  ha-e 


(54)  &k(*;  *,+)*  ) 


0<x<i' 


for  X  >  i.  , 


where 


fa  %ji  -  p-o/  f/nM'*  >]/„.„■ 


is  the  residue  of  the  integrand  of  (43)  at  the  K  th-order  pole  at 
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Clevly  r(s+fi+j+/)  * 

>"*j-0(A+'V 

where  again  an  'empty'  product  should  be  interpreted  as  unity.  Hence 


r(p+sj0>4"jj*  (h*t+to(b**vmt)"'Lb+t]  "  ty+w+ij  JUb+t+kJ  t 

nr 


(56)  X*te*JJ  •  ctryfe f rt" r,lw*oJ[ ^ 


r* 


Denoting  by  arjOj  the  function  in  brackets  in  (56)  we  obtain: 


(57)  (rJ(tj  =  Lbritkj'  yr^fitl+lj 

•  A-t  ' 


-K/r,V 


and 


log  £■;(*)  =*  _  xiorx  -log*  +  h 

'  A*® 


The  successive  logarithmic  derivatives  of  lO  are  given  by 

4///J  *"'ogx  -  J-'*  -*.+if>4-cf4S)  -  Z.CP+l.itj'' 

(50]  J  J  J  *■  A*e  /•< 

/‘"=  L-5"Wr""'<  «/>.>/;  -  ^‘(nxtfSj <■ e>rr!Z.Lh*» 

4  kt* 

and  since  (55)  are  to  be  taken  at  j  *  -/.  we  have 

^•r(  vV  t*v)  ~L-o  Ky  *  '/  r  \  i S  J  /  '  , 

(50)  4;c-/»v/  =  -l-Jt:*  * "■Y'Li’jJ  tK£t, *  :* 

-  <rOrr'r\^j-irt,U^r^  -nfci-jj+nr.'  Li~LtJ 


*"■»») 


-)/ 

and  finally 
(60) 


^  *  ■  •  tr  //  ’  w  #  '  *  f  " £ -y  y 

.  rw*;  j.  y  xr/^rJ -L~  7 

I  rHn)  LA',J'pZ  AV  J 


*  * 


where  the  arguments  of  -function  are  given  by  (5S). 


The  convergence  of  the  infinite  series  (60)  follows  from  a  similar 


argument  as  that  applied  in  the  case  of  the  infinite  series  (1*0). 
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4,  Distribution  of  products  of  Independent  random  variables  having 
the  same  Beta  distributions  with  different  scale  parameters. 

The  generalization  of  the  above  results  to  the  case  important  In 
practical  applications  vhen  y  ■  -X*.  and  the  p.d.f.  of  the 

i  th  factor  has  the  p.d.f. 

*4-  p  C*<x''  r3'ty) 

is  Immediate.  The  p.d.f  fa(yj  of  such  a  product  Is  equal  to.’ 

(6i)  *  */  /J*  (*!*».•  •- z*y J  ^  ) 

and,  consequently  the  probability  distribution  function  Is  equal  to.' 

8* ft,*]. 

Indeed^  since 

ve  have  . * 

•  -w''*'  n  1*k( tr.uhj' 

*  *,**••«*  ^l[ hL^^  -^yy/i,^/ i ] 

r  M  [  *,<*>.  **  1 J  j 

which  completes  the  proof  of  (6l)  In  view  of  the  uniqueness  property  of  the 
Mellin  transform  (cf.  (3),  p.  104,  Th.  2). 


r 


v 


*  I 


(1A) 


(°A) 


(JA) 
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APPSTDIX 

(a)  The  evaluation  of  the  Integral  (17)  by  contour  integration  leading 
to  formula  (10)  la  achieved  in  the  following  way: 

We  define  two  cloae  contours  and  . 

o C4.  is  composed  of  the  vertical  segaent  of  the  integration  path  of  (17) 
between  C-i4  and  c  +  i'4-  and  of  a  half -circle  %A.ot  radius  A-  lying 
to  the  right  of  this  segment  and  having  this  segment  aa  its  diameter. 
Similarly 

4  la  composed  of  the  same  vertical  segment  and  a  half -circle  of 

radius  A  lying  to  the  left  of  this  segment  and  having  this  segment  as 


-  >v 
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is  analytic  on  the  contours  Cfi  and  and  its  interiors  with  the 

exception  of  ^  poles  inside  the  contour  we  have 

f  fit)**  -  0  f  ft*)**  *  Z 

■4f  ^4  J~* 

where  £(x;  njj j  iB  the  residue  of  f(ij  at  the  pole  S  -  -  /  (>0,1,...  f-t  ). 

It  remains  only  to  prove  that  the  integral  along  the  half -circle  ^4  tends 

to  zero  when  A  tends  to  infinity  in  the  case  of  X,  >  1  and  that  the  integral 

along  the  half -circle  5^,  tends  to  zero  when  A  tends  to  infinity  in 

the  case  of  0  <  x  <  1. 

(b)  Let  X  ^  1.  Then  y, 

But  in  view  of  cos  being  positive  we  have 

J  C4  k.  +  Acnf  t  A 1  I  *  /  (e+4Jh+  lAlc+4)tnf+41'l  ?  A 

tt‘  a-*V 

and  the  product  //^  in  (5A)  is  dominated  by  A  .  Hence 

/  /  fUXi  /  *  2  A"**  fit"*  X~  A  A  $ 

*4  9  w 

Since  cos  f  >  0  for  0  <  /z  and  X  >,  I  J  i  —"2  . 

In  view  of  (6A)  the  integral  along  the  half -circle  tends  to  zero  when  A  tends 
to  infinity  thus  completing  the  proof  of  the  Becond  line  of  (16). 
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(c)  Let  0  <  x  <1.  Then 

r  ,%rtX_Cc+A<»t+AiJi»+)V  \  /a 

/  «  J  X  7T (C-ti+Arnfi-AJn+J  Alt  V* 

*‘A  % 

But  ^ 

jc±k  +  Arlfl  f>  +  Ai  Jin  +I=(CC  +kj*+  %A( Hkjcai  f  + A  J 

end  this  Is  greater  than 

/(c*Ajl-  2A(C+*)+A1/*  *  I  A- L-k!  >  4-C-f+f  >A-f 

since  C  >  0  could  be  made  less  than  unity,  and  It  should  be  borne  In 

mind  that  we  have  assumed  A*>  <y ,  Hence  the  product  In  (7A)  Is  dominated 
(•  r** 

by  the  (tv  and  r 

Hk 

*  ***€*+ 

The  Integral  above  can  be  written  {x^+dt  and  since  cos  A  >0 
and  0  <  X  *$:  1  /s^ows  that  the  left-hand  side  of  (8A)  tends  to  zero 
when  4-  tends  to  Infinity,  thus  completing  the  proof  of  the  first  line  of  (16). 


(d)  The  evaluation  of  the  Integral  (3*0  by  contour  Integration  leading 
to  fonnula  (3^)  proceeds  in  a  similar  way.  Notice  first  that  the  Integrand 
of  (3U) 

[*”  1 K 

fct)  =  x  Lri**-+)J 

is  analytic  with  the  exception  of  the  Infinity  of  poles  for  5  ■  0,-1, -2.... 
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Let  X  1,  Then  on  the  contour  as  defined  above  by  (lA)  ve  have 
S-  U.  +  «>  with  U,  >  0  and 

l  r(t)  I  ICiLtlJ  I  m  /  rf/y-~  / 

I  rjuy)l  I  r(f*f)  I  • 


(ioa) 


But  for  complex  X(see  (1], formula  6.1.25)  and  X  ^  0,-l,-2.... 

jr(jjl=  r(*)  ft  ['  ] 


(UA) 


so  that  00 

/  1  _  7/ 

/  r(Sr^j  I  ~  r(ut$J  r*o 


/  +Cu%*r)i" 

- -  -y-r  -  - 

/  ’+(14.41 4  r)*- 


i 


since  with  ^  >  1  all  the  factors  of  the  Infinite  product  above  are  smaller 
than  unity;  hence 


by  the  same  argument  as  that  applied  to  the  evaluation  of  (5A)  and  (6A) 
we  see  that,  for  x  ^  1,  the  integral  of  (9A)  along  the  half -circle 
tends  to  zero  when  A  tends  to  infinity  thus  completing  the  proof  of  the 
second  line  of  (35). 
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(ldA) 


(e)  To  prove  the  first  line  of  (35)  (for  0  <Jt  <  1)  ve  define  a 
contour  slightly  different  from  that  defined  by  (2A);  it  is  shown 

in  Pig.  3.  Let 

be  composed  of  the  vertical  segment  of 
the  integration  path  of  (3M  between 
and  c *  A- 1  f  of  two  horizontal  segments  Joining 
in  the  5  -plane  points  -AiAi  with  t  +  Ai 
and  -  A  -  At  with  C-  At  and  finally  of  the 

vertical  segment  Joining  point  -  A  -  At  with  Figure  3 

-A  *■  A’i  ,  To  avoid  the  poles  of  function  (9A)  we  chose  a  sequence  of 
/f- values  as  Am  i  and  we  are  going  to  show  that  the  integral  of 

(9A)  along  the  path  composed  of  segments  marked  (i)  -  (v)  in  Fig,  3 
tends  to  zero  when  A  tends  to  infinity  assuming  the  values  A~  m  yn+\  , 


-  A (« i 

liic  ♦  A  • 

• 

1 

• 

* 

1 

1 

1 

1 

1 

< 

« 

f 

1  ( 

i 

1 

t 

i 

1 

« 

i 

9 

• 

1 

1 

| 

♦ 

m  m  m  m  i 

-  A-Ai  l /*i 

V) 

(f)  In  the  integral  along  (v)  the  real  part  of  J  «  U.-AI  is 

positive  and  by  (12A)  the  integral  is  dominated  by 

I  I  I  s  '  fC  ~t  I  /  /"*““■  Am,  -  L  .  * 

/  J  K  1  *■  fcXr  i*  f*  -e.,r 

'  I'J  10 

which  tends  to  zero  when  A-  tends  to  infinity;  the  same  applies  to  the 
integral  along  (l). 


(g)  To  show  that  the  Integrals  along  segments  (ii),  ( 111 (iv) f  where 
the  real  part  of  S  is  negative,  have  the  same  property  notice  that  for  any 
non-integer  z  (cf.  for  example  (lj,  formula  5.1.17); 
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(HA) 


r (zj  r(t  -  zj  *  j(tx(n.)  / 


so  that 


r(-  2  +  /  J  =  r(z)  JtK(nj 


and 


IT 


r(~z  +  f)  -  f/T^77) 


Writing  in  the  above  —  S  for  z.  ve  have 

r(s+0  _  r(-s-}+0 


But  in  the  Integrals  which  we  are  going  to  evaluate  the  real  part  of  -J 
is  positive  and  we  can  apply  (10A)  which,  in  view  of  f  >  1  yields 


(15A) 


/  rfu/i  i  j 


(h)  In  the  Integral  along  (ill) 

,  ,  ,f  »*»• 

I  ffiijM  <IJX  rV4.*«.v/  (-*"") 

(Jit)  “A 

It  is  known  that  if  z  »  k^/V  then 


J  ■  and  in  view  of  (15A) 

7<r  XV  /  *-•*!]  I  (A 


i 


!)!■»*/  m  [Jin\  +  (cf  # 


for  example  [l],  formula  4.3.59) 
<• 


so  that 


4 
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Recalling  that  <4  ■  *H^X  we  have  //*i  K k  ■  i  and  the  Integrand  above 
la  not  greater  than  unity;  the  Integral  la  dominated  by  which 

tends  to  zero  when  A  tends  to  Infinity. 


(l)  In  the  Integral  along  (iv)  S  ■  U-t  A  and 


/  //'Mi 


1  x  ry *-**+?)  fc-tA)  . 


-A 


In  view  of  (15A)  this  is  dominated  by 

,  [caVtoli 

'  u 

and  this  tends  to  zero  when  ^  tends  to  infinity;  the  same  applies  to 
the  integral  along  (ii). 

(j)  When  A  m  Z  the  contour  includes  in  its  interior  *n,+  1  poles 
of  function  (9A)  and  with  r\  tending  to  infinity  the  integral  (34)  is  equal  to 

r\i+o  6  m  ZL  Kc*,  *-jj) , 

9H  -?O0  J  ~0 


which  is  the  infinite  convergent  series  (40). 

(k)  The  evaluation  of  integral  (44)  by  contour  integration  proceeds 
in  a  similar  way  as  that  of  integral  (17)  and  of  integral  (43)  (in  the  case 
of  not  being  an  integer)  in  a  similar  way  as  that  applied  in  the 

evaluation  of  integral  (34).  In  fact  the  task  is  easier  with  the  additional 
factor  J  appearing  in  the  relevant  integrals. 


Summary 


The  Important  problem  of  finding  the  probability  density  function  of 
the  product  of  a  number  of  independent  identically  distributed  random 
variables  was  Investigated  by  Springer  and  Thompson  [8],  [9]  who,  by  the 
application  of  the  Mellin  transform,  obtained  the  ansver  In  the  case  of 
the  normal  and  the  Cauchy  random  variables.  Further  discussion  by  the 
author  of  this  report  (7]  has  dravn  attention  to  the  fact  that  the  results 
obtained  by  them  are  easily  applicable  to  the  cases  of  exponential,  gamma 
and  Welbull  distributions  and  that  the  probability  distribution  functions 
can  be  also  obtained  directly  by  the  application  of  the  Mellin  transform 
without  integrating  the  relevant  formulae  for  the  probability  density 
functions. 

In  the  present  paper  it  is  shown  how  these  methods  can  be  applied 
to  the  study  of  beta  random  variables;  this  is  of  practical  Importance 
since  most  physical  devices  with  which  an  engineer  is  dealing  have  finite 
characteristics  and  the  results  discussed  in  [8],  (9],  [7l  were  mainly 
concerned  with  random  variables  having  infinite  ranges. 
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The  problem  of  finding  the  probability  distribution  of  a  product 
of  a  number  of  identically  distributed,  Independent  random  variables  had  been 
solved  by  an  application  of  the  Me.llin  transform  for  normal  and  Cauchy  distribu¬ 
tion'  (Springer  and  Thompson)  and  for  exponential,  Gamma  and  Welbull  distributions 
(Lomnicki).  The  present  paper  shovs  that  it  can  be  solved  by  '•imllar  methods  for 
B'f.  distributions.  This  Is  of  rracticel  importance,  inre  m.06t  physical  quantitle. 
with  which  an  engineer  is  dealing  have  finite  range  ,  while  all  the  distributions 
previously  studied  had  infinite  ranges. 
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